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We show that a simple first-principles correction based on the difference between the singlet- 
triplet CO excitation energy values obtained by DFT and high-level quantum chemistry methods 
yields accurate CO adsorption properties on a variety of metal surfaces. We demonstrate a linear 
relationship between the CO adsorption energy and the CO singlet-triplet splitting, similar to the 
linear dependence of CO adsorption energy on the energy of the CO 2n* orbital found recently 
[Kresse et ai, Physical Review B 68, 073401 (2003)]. Converged DFT calculations underestimate the 
CO singlet-triplet excitation energy AEs-t, whereas coupled-cluster and CI calculations reproduce 
the experimental AEs-t- The dependence of -Echom on A_Es-t is used to extrapolate J5 c hcm for the 
top, bridge and hollow sites for the (100) and (111) surfaces of Pt, Rh, Pd and Cu to the values that 
correspond to the coupled-cluster and CI Ai?g_T value. The correction reproduces experimental 
adsorption site preference for all cases and obtains i5 c hcm in excellent agreement with experimental 
results. 



I. INTRODUCTION 

The chemisorption of carbon monoxide on transition 
metal surfaces is regarded as a prototypical system for 
the study of molecule-surface interactions and has been 
intensely studied theoretically and experimentally 0. 
For first-principles theoretical studies, density functional 
theory (DFT) with the generalized gradient approxima- 
tion (GGA) to the exchange-correlation functional has 
emerged as the method of choice, and DFT-GGA stud- 
ies have significantly advanced understanding of surface 
phenomena 

Despite these successes, two of the most basic proper- 
ties of CO-metal surface interactions - the chemisorp- 
tion energy and the preferred adsorption site, cannot 
be reliably predicted by DFT calculations. Theoretical 
CO chemisorption energies obtained by the widely used 
PW91 and PBE functionals are significantly higher 
than experimental values, sometimes by as much as 0.4 
eV @ (30%). The RPBE functional [| does improve the 
adsorption energetics, but at the expense of lower accu- 
racy in metal lattice constants and a severe underestima- 
tion of surface energies. Even more importantly, neither 
PBE nor RPBE calculations can correctly predict the rel- 
ative energetics of the high symmetry sites, favoring the 
more coordinated bridge and hollow sites over the top 
site, and resulting in a wrong site preference in a num- 
ber of systems. This was first discussed by Feibclman 
et al. 3 where a variety of DFT-GGA methodologies 
and codes predicted adsorption at the fee or hep hollow 
site to be preferred over the experimentally preferred top 
site adsorption on the Pt(lll) surface at low coverage. 
Since then this "puzzle" has been addressed y| 01 UH LU| 
and the DFT-GGA inaccuracy traced to the incorrect de- 
scription of CO electronic structure and bond-breaking. 

Grinberg et al. showed that the inaccuracy in site pref- 
erence was due to the poor treatment of CO bond break- 
ing by GGA functionals. Analyzing the CO electronic 



structure, Gil et al. found that PBE and B3LYP [T1[T^| 
functionals place the unfilled CO 2tt* orbital too low in 
energy. This makes it too close in energy to the metal ri- 
band, which results in an unrealistic strengthening of the 
2ir*- d-band bonding interaction. More recently, Kresse 
and co-workers [Tj used DFT calculations of varying ac- 
curacy and showed that a linear relationship exists be- 
tween the difference of top and hollow site chemisorption 
energies for CO on Pt(lll) and the gas-phase energy of 
the CO 2tt* orbital. 

By using a GGA + U type functional, Kresse et al. \H\ 
were able to adjust the energetic position of the gas-phase 
CO 2ir* orbital, restoring the correct prediction that top 
site adsorption is preferred on Pt(lll). However, the 
"correct" value of U is not known a priori. Further- 
more, they studied adsorption on Pt(lll) and did not 
address CO on other metal surfaces. In this paper, we 
demonstrate that a linear relationship exists between the 
CO chemisorption energy and the CO singlet-triplet ex- 
citation energy for the top, bridge and hollow sites on 
a variety of metal surfaces. Unlike the energetic posi- 
tion of the 2n* orbital, the CO singlet-triplet excitation 
energy is rigorously well-defined and is accurately com- 
puted by coupled-cluster |14J and configuration interac- 
tion (CI) quantum chemical calculations ^5|. Extrapola- 
tion of -Echcm values to the correct CO singlet-triplet ex- 
citation energy relies only on first-principles calculations, 
and yields chemisorption and site-preference energies in 
excellent agreement with experiment for all systems stud- 
ied. 



II. METHODOLOGY 

Calculations are carried out using the PBE GGA 
exchange-correlation functional |5| and norm-conserving 
optimized pseudopotentials 01 with the designed nonlo- 
cal method |l7| for metals. All pseudopotentials were 
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constructed using the OPIUM pseudopotential pack- 
age . CO chemisorption is modeled at 1 /4 monolayer 
coverage on five layer slabs, separated by vacuum, with 
allowed relaxation in the top two layers. All calculations 
are done, and tested to be converged, using a 4 x 4 x 1 
grid of Monkhorst-Pack fc-points 19] . The Kohn-Sham 
orbitals are expanded in a plane- wave basis set truncated 
at either 81 or 50 Ry, the higher cutoff being required in 
calculations using C and O pseudopotentials with small 
real-space cutoffs. We calculate the chemisorption en- 
ergy for CO adsorbed on the (111) and (100) surfaces of 
Pt, Rh, Pd, and Cu, in three high symmetry sites: top, 
bridge, and hollow. On the (111) surfaces, adsorption at 
the fee hollow site is not reported, since it is well known 
that differences between calculated fee and hep hollow 
site adsorption are negligible |2(l. l2l|. 

We repeat these calculations using different sets of C 
and O pseudopotentials described in Table I. We use 
three sets of C and O pseudopotentials to calculate 
E c hem at each site on each surface. To evaluate whether 
three data points are adequate to describe trends in the 
chemisorption energies, we expand the number of C and 
O pseudopotential sets to five, and repeat the calcula- 
tions for the hollow site on Pt(lll). The slope for 
versus AEs-t is unchanged on going from three points 
to five, as is the goodness of fit. Based on this, the re- 
mainder of our results use the first three pseudopotential 
sets. 



III. 



RESULTS 



For CO on Pt(lll), our calculations show a linear re- 
lationship between E c ^ m and CO Ai?s_T excitation en- 
ergy (Figure 1), similar to the linear relationship between 
CO/Pt(lll) Ai? top -fcc site preference energy and the en- 
ergy of the CO 2ir* for the Pt(lll) adsorption observed 
by Kresse and coworkers A linear fit also describes 
well the behavior of chemisorption energies on the seven 
other substrates included in the present work. Since the 
CO triplet state is produced by an excitation of an elec- 
tron from the 5cr orbital to the 27r* orbital, the singlet- 
triplet excitation energy is closely related to the 5a-2ir* 
gap and to the position of the 2ir* orbital. This gives rise 
to a similar linear dependence of -B^hem on the energy of 
CO 2?r* orbital and on AE's-t- 

Coupled-cluster [14( and CI [l5( quantum chemical cal- 
culations accurately reproduce the experimental AE^-t 
of 6.095 eV |2^. On the other hand, regardless of 
the pseudopotential set, our DFT-GGA calculations al- 
ways predict a Ai?g_T that is too small. The correct 
chemisorption energy E^" m can be obtained by using 
the relationship between -E c hcm and Ai?s-T and extrap- 
olating to the CI Ai?s-T value: 
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where AE$]_ T and AE^St are respectively the CI and 



GGA CO singlet-triplet excitation energies, and SE^^ 
I SAE^t is the slope of the fit of -E c hem versus AE$-t- 

A universal feature of the corrected chemisorption 
energies in Table II is that they all indicate weaker 
chemisorption than the corresponding uncorrected val- 
ues, with the ££hcm values for the preferred site demon- 
strating much better agreement with experimentally de- 
termined adsorption energies. For example, the £"™^ for 
the experimentally seen CO/Pd(lll) hollow site is 1.96 
eV, as compared to 1.47-1.53 eV obtained by temperature 
programmed desorption (TPD) measurements 23] . This 
rather large 0.46 eV error is eliminated by the use of the 
extrapolation, with ££hcm of 1.60 eV, in very close agree- 
ment with experimental results. For the experimentally 
observed top site on Cu(lll) the results are less dramatic 
but still noticeable, with E%^£ of 0.746 eV changed to 
^chorn °f 0-621 eV in better agreement with the experi- 
mental value of 0.49 eV [24| . Overall, comparison of PBE 
and corrected results shows an improvement from 0.38 eV 
(30%) average overestimation for the PBE functional to 
0.16 eV (13%) average error for our corrected results. 

Examination of the data in Table II shows that there 
is a strong correlation between the magnitude of the cor- 
rection and the chemisorption site. This is due to the 
different strengths of the metal-CO interactions in differ- 
ent local geometries. While the chemisorption bond is 
formed through both a donation and 7r* back-donation, 
the contribution of 7r* back-donation to the adsorption 
bonding in the systems considered here is dominant be- 
cause the fillings of the late transition metals studied here 
are greater than half [23 ■ The back-donation mechanism 
is strongly enhanced by going f rom top site to polycoor- 
dinated adsorption sites [OjESIIEIIW- Accordingly, the 
incorrect DFT-GGA AEs~t (or the incorrect placement 
of the 2ir* orbital) will require the smallest correction 
for £7 c hem °f the top site, followed by the bridge site, the 
three-fold hep hollow site on the (111) surfaces and the 
four-fold hollow on the (100) surfaces. This ranking is ev- 
ident in the tabulated results. In the chemical language, 
the CO bond weakening is smallest for the top site and 
largest for the four-fold hollow site. 

The DFT-GGA errors in the prediction of the preferred 
site are a direct outcome of the unequal treatment of 
the CO bond weakening at the top and hollow sites due 
to the unrealistically small Ai?g_x and low 2ir* energy. 
Using our first-principles extrapolation to eliminate the 
CO bond weakening errors, our corrected DFT results 
give the highest value of E^"^ for the experimentally 
observed sites in all cases. For Pt(lll), our raw DFT 
data show an incorrect site preference with an energy of 
0.076 eV. The corrected e nerg ies agree with experimen- 
tal site preference [2H, I29I l30j] with an energy difference 
of 0.163 eV. Likewise, our raw DFT data for Rh(lll) is 
disparate with experimental site preference 0, llj while 
our corrected energies are in agreement. For Cu(lll) the 
use of the correction gives the experimental site prefer- 
ence 0, H3 ! though the corrected DFT results predict 
small (< 0.1 eV) differences between the top, bridge and 
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hollow sites. Both the raw DFT and the corrected re- 
sults for Pd(ll l) agree with the site preference observed 
by experiment jUI^ISl- Our corrected results for the 
(100) surfaces agree with experimentally observed pre- 
ferred adsorption sites for Pt [3^, |33, El] , Rh [3^, > 
Pd [alliil, and Cu jHEI. 

A less exact but simpler correction can be extracted 
from our data and applied to any CO /metal surface sys- 
tem. The data in Table II show that for any given site 
the SE^m I SAEgSr values and the consequent cor- 
rection energy are fairly constant with scatter of about 
0.1 eV across a range of systems. 

As discussed above, the similarities for the same site 
and the differences among the sites are consequences of 
the different degrees of CO bond weakening. The degree 
of CO bond weakening can be estimated from the fre- 
quency shift of adsorbed CO relative to the gas phase CO 
molecule, which can be easily calculated for any system. 
Then the frequency can be compared to those of top (typ- 
ically 2000-2100 cm- 1 ), bridge (1850-1950 cm" 1 ), hep 
(1750-1800 cm" 1 ) or four-fold hollow (1600-1700 cm- 1 ) 
sites and a corresponding correction applied. 

For example, for CO adsorption on Ni(lll), DFT- 
GGA calculations with PBE or PW91 functional find 
the preferred site to be hep or fee [44j , in agreement 
with experimental results [4fij . However, the E c hem are 
in 1.9-2.0 eV range, in contrast to experimental £^hem of 
1.12-1.55 eV. Since the hollow site CO/Ni(lll) stretch 
frequency of 1800 cm -1 |5| is similar to the 1830 cm -1 
frequency for CO on Pd(lll) hollow site [2l|. we expect 
the -Echcm error to be similar, and the corrected value of 
-Echcm is about 1.55 eV, in much better agreement with 
experimental results. A graphical representation of the 
suggested correction to E c h cm as a function of CO stretch 
frequency is given in Figure 2. 



TABLE I: Pseudopotential details. Core radii are in a , plane 
wave cutoffs q c in Ry. All PSPs were created from the s 2 p 2 
reference configuration for carbon and the s 2 p 4 reference con- 
figuration for oxygen. For each pseudopotential set, results 
from gas-phase molecule calculations for the 2n* energy, as 
well as the 5a-2-K* gap and the singlet-triplet energy are given, 
all in eV. 
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PSP 1 0.94,1.09 


81,81 


-2.10 


6.91 


5.35 


PSP 2 1.60,1.49 


47,50 


-1.90 


7.01 


5.53 


PSP 3 1.70,1.49 


30,50 


-1.61 


7.35 


5.84 


PSP 4 1.65,1.49 


42,50 


-1.94 


7.04 


5.56 


PSP 5 1.70,1.49 


39,50 


-1.87 


7.09 


5.61 
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IV. CONCLUSION 



We have shown that the chemisorption energies of CO 
adsorbed on metal surfaces depend linearly on gas-phase 
CO singlet-triplet splitting. The difference between the 
high level quantum chemistry coupled-cluster/CI and 
DFT-GGA singlet-triplet excitation energies can then be 
used to extrapolate to chemisorption energies with the 
CO error removed. The corrected values are in good 
agreement with experimental results. The correction also 
eliminates the GGA errors in site preference. We find 
a strong correlation between the amount of CO bond 
breaking and the correction magnitude. This suggests 
that an estimate of the GGA error due to the incorrect 
description of CO electronic structure can be readily ob- 
tained through the frequency shift of adsorbed CO at 
any site on a metal surface. The demonstrated method 
should be applicable to various adsorption systems where 
the charge transfer responsible for chemisorption is sen- 
sitive to the adsorbate electronic structure. 
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TABLE II: Results of linear regression of chemisorption en- 
ergy versus singlet-triplet splitting energy. The correlation co- 
efficient for each fit is in parenthesis next to each slope. The 
DFT/GGA values for the chemisorption energies are given, 
along with the corrected energies obtained by extrapolation. 
For convenience, the value of the difference between top and 
hollow site adsorption energies, E-t—h, an d how it evolves with 
the magnitude of correction, A, are also listed. Positive value 
of E t -h indicates that top site is preferred. For each substrate, 
the site found to be preferred by experiment is marked with an 
asterisk (*). Experimental values for E c h em are given in paren- 
theses next to E^gx for Pt(lll) [pJil Ej , Rh(lll)JHp3, 

Pd(in) J2I J3a, Cu(iii)jMT2arPt(ioo) |M ETUI, 

Rh(100) mEQl7Pd(100) ilQ, and Cu(100) glEl 
Site Slope £ c G hem £c c h °c r m A 



Pt(lll) 


top* 
bridge 
hep 
Et-h 


-0.211(0.989) 
-0.435(0.999) 
-u.ooz(u.yyy j 
0.321 


Rh(lll) 


top* 
bridge 
hep 
Et-h 


-0.259(0.975) 
-0.456(0.997) 
-u.ooy(u.yyD j 
0.300 


Pd(lll) 


top 
bridge 
hep* 
Et-h 


-0.185(0.977) 
-0.399(0.977) 

/o nnt;^ 

-u.ooo (u.yyo, 
0.350 


Cu(lll) 


top* 
bridge 
hep 
Et-h 


-0.169(0.981) 
-0.329(0.951) 
-0.375(0.993) 
0.206 


Pt(100) 


top* 
bridge* 
hollow 

Et-h 


-0.212(0.975) 
-0.422(0.996) 
-0.607(0.999) 
0.395 


Rh(100) 


top* 
bridge* 
hollow 

Et-h 


-0.246(0.989) 
-0.427(0.999) 
-0.651(0.999) 
0.405 


Pd(100) 


top 
bridge* 
hollow 
E t -h 


-0.196(0.980) 
-0.384(0.995) 
-0.583(0.999) 
0.387 


Cu(100) 


top* 
bridge 
hollow 

Et-h 


-0.170(0.996) 
-0.286(0.991) 
-0.523(0.999) 
0.353 



1.717 1.560(1.43-1.71) -0.157 

1.758 1.433 -0.325 

1.793 1.397 -0.396 

-0.076 0.164 0.240 

1.866 1.673(1.43-1.65) -0.193 

1.920 1.581 -0.339 

2.059 1.644 -0.415 

-0.193 0.030 0.223 

1.385 1.247 -0.138 

1.784 1.487 -0.297 

1.962 1.602(1.47-1.54) -0.360 

-0.577 -0.355 0.222 

0.746 0.621(0.46-0.52) -0.125 

0.822 0.576 -0.246 

0.889 0.610 -0.279 

-0.143 0.011 0.154 

1.954 1.796(1.62-2.18) -0.158 

2.139 1.824(1.62-2.18) -.315 

1.698 1.246 -0.452 

0.256 0.551 0.295 

1.905 1.723(1.24-1.65) -0.182 

2.092 1.774 -0.318 

2.087 1.603 -0.484 

-0.182 0.120 0.302 

1.494 1.348 -0.146 

1.927 1.641(1.3-1.71) -0.286 

1.937 1.503 -0.434 

-0.443 -0.155 0.288 

0.830 0.703(0.55-0.57) -0.147 

0.834 0.620 -0.214 

0.831 0.441 -0.477 

-0.001 0.262 0.263 
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FIG. 1: £ c h h c c p m for CO on Pt(lll) versus AE S -t of CO for 
five pseudopotential sets (circles) and as determined by CI 
calculations (square). Similar clear linear relationships are 
also obtained for other sites and metal surfaces. For both the 
most accurate pseudopotential set as well as for the CI data 
point, dashed lines guide the eye to the corresponding val- 



ues of EJ^ m . First-principles extrapolation procedure: E^^ 



and AEg3r~ values are plotted and fit to a line (solid). This 
line is extended to the abscissa representing AE^Lt- The cor- 
responding ordinate value is AE^ll^. Chemisorption on the 
hep site of Pt(lll) is used in this example. 




1600 



1700 



1800 



1900 



2000 



2100 



-L 



CO Stretch Frequency (cm ) 



FIG. 2: Graphical guide for estimating E^Sm from CO stretch 
frequency. Correction energy as a function of the adsorbed 
CO stretch frequency provides a simpler, more user-friendly 
method of applying the correction. Lines represent bound- 
aries for suggested correction values, average correction val- 
ues for high-symmetry sites for our highest quality DFT-GGA 
results in this study are shown as circles. 
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